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Genetically, the classical evolution of the inflaton has a brief fast roll stage that precedes the 
slow roll regime. The fast roll stage leads to a purely attractive potential in the wave equations of 
curvature and tensor perturbations (while the potential is purely repulsive in the slow roll stage). 
This attractive potential leads to a depression of the CMB quadrupole moment for the curvature 
and B-mode angular power spectra. A single new parameter emerges in this way in the early 
universe model: the comoving wave number fci characteristic scale of this attractive potential. This 
mode fci happens to exit the horizon precisely at the transition from the fast-roll to the slow-roll 
stage. The fast-roll stage dynamically modifies the initial power spectrum by a transfer function 
D(k). We compute D(k) by solving the inflaton evolution equations. D(k) effectively suppresses 
the primordial power for k < ki and posseses the scaling property D(k) — 9(k/ki) where 't'(x) is 
an universal function. We perform a MCMC analysis of the WMAP and SDSS data including the 
fast-roll stage and find the value fci = 0.266 Gpc -1 . The quadrupole mode kQ — 0.242 Gpc -1 exits 
the horizon earlier than fci, about one-tenth of an efold before the end of fast-roll. We compare 
the fast-roll fit with a fit without fast roll but including a sharp lower cutoff on the primordial 
power. Fast-roll provides a slightly better fit than a sharp cutoff for the TT, TE and EE modes. 
Moreover, our fits provide non-zero lower bounds for r, while the values of the other cosmological 
parameters are essentially those of the pure ACDM model. We display the real space two point 
C TT (8) correlator. The fact that kQ exits the horizon before the slow-roll stage implies an upper 
bound in the total number of efolds N to t during inflation. Combining this with estimates during 
the radiation dominated era we obtain Ntot ~ 66, with the bounds 62 < N to t < 82. We repeated 
the same analysis with the WMAP-5, ACBAR-2007 and SDSS data confirming the overall picture. 
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I. INTRODUCTION AND RESULTS 



The Standard (Concordance) Model of the Universe explains today a wide set of cosmological and astronomical 
measurements performed over a large variety of wave-lenghts and observation tools: large and small angular scale CMB 
observations, light elements abundances, large scale structure observations (LSS) and properties of galaxy clusters, 
Hubble Space Telescope measurements on the Hubble constant, supernova luminosity/distance relations (acceleration 
of the today universe expansion), and other measurements. The concordance of these data imply that our universe is 
spatially flat, with gravity and cosmological perturbations described by Einstein General Relativity theory. WMAP 
data give a strong support to the Standard Model of the Universe. 

Inflation was introduced to solve several outstanding problems of the standard Big Bang model [1[ and has now 
become an important part of the Standard Model of the Universe. At the same time, it provides a natural mechanism 
for the generation of scalar density fluctuations that seed large scale structure, thus explaining the origin of the tem- 
perature anisotropics in the cosmic microwave background (CMB), as well as that of tensor perturbations (primordial 
gravitational waves) @, [l6| . 

The horizon and flatness problems are solved provided the universe expands for more than 62 efolds during inflation. 
This is achieved within slow-roll inflation where the inflaton potential is fairly flat. 

Although there are no statistically significant departures from the slow roll inflationary scenario at small angular 
scales (I > 100), the WMAP data again confirms the surprinsingly low quadrupoles Cj T and C% E [l]-[9[ and suggests 
that it cannot be completely explained by galactic foreground contamination. The low value of the quadrupole has 
been an intriguing feature on large angular scales since first observed by COBE/DMR [ij], and confirmed by the 
WMAP data @]-[|. 

In order to asses the relevance of the observed quadrupole suppression in the ACDM model, we determine in the best 
fit ACDM model the probability to observe the quadrupole 20% below the theoretical mean value. This probability 
turns out to be only ~ 0.06. This small probability supports the necessity for a cosmological explanation of the 
quadrupole depression beyond the ACDM model. 

Generically, the classical evolution of the inflaton has a brief fast roll stage that precedes the slow roll regime. 
The fast roll stage leads to a purely attractive potential in the wave equations of curvature and tensor perturbations. 
Such potential is a generic feature of this brief fast roll stage that merges smoothly with slow roll inflation. This stage 
is a consequence of generic initial conditions for the classical inflaton dynamics in which the kinetic and potential 
energy of the inflaton are of the same order, namely, the energy scale of slow roll inflation. During the early fast roll 
stage the inflaton evolves rapidly during a brief period, but slows down by the cosmological expansion settling in the 
slow roll stage in which the kinetic energy of the inflaton is much smaller than its potential energy. 

As shown in ref. @, @] the attractive potential in the wave equations of curvature and tensor perturbations during 
the fast-roll stage leads to a suppression of the quadrupole moment for CMB and B-mode angular power spectra. 
Both scalar and tensor low multipoles are suppressed. However, the potential for tensor perturbations is about an 
order of magnitude smaller than the one for scalar fluctuations and hence the suppression of low t tensor perturbations 
is much less significative [1, 0] ■ 

The observation of a low quadrupole [H H, and the surprising alignment of quadrupole and octupole [l(| El 
sparked many different proposals for their explanation [l2T |. 

The fast-roll explanation of the quadrupole does not require to introduce new physics neither modifications of the 
slow- roll inflationary models. The only new feature is that the quadrupole mode should exit the horizon during the 
generic fast-roll stage that precedes slow-roll inflation. 

A single new parameter emerges dynamically due to the fast-roll stage: the comoving wave number k%, characteristic 
scale of the attractive potential felt by the fluctuations during fast-roll. The fast-roll stage modifies the initial power 
spectrum by a transfer function D{k) that we compute solving the classical inflaton evolution equations [see fig. 3]. 
D{k) effectively suppresses the primordial power for k < k\ and posseses the scaling property D{k) — ^(fc/fci) where 
^(x) is an universal function. D(k) has a main peak around ku — 1-9 k\ and oscillates around zero with decreasing 
amplitude as a function of k for k > ku- D(k) vanishes asymptotically for large k, as expected. 

We report here the results of a MCMC analysis of the WMAP-3, small-scale CMB and SDSS data including the 
fast-roll stage and find the value k\ — 0.266 Gpc -1 . This mode k\ happens to exit the horizon precisely at the 
transition from the fast-roll to the slow-roll stage. The quadrupole mode kq — 0.242 Gpc -1 exits the horizon during 
the fast-roll stage approximately 1/10 of an efold earlier than k\. We compare the fast- roll fit with a fit without 



3 



fast roll but including a sharp lower cutoff on the primordial power. Fast-roll provides a slightly better fit than a 
sharp cutoff for the Cj T , Cj E and Cf E coefficients. Besides reproducing the quadrupole depression, the fast roll fit 
accounts for the oscillations of the lower multipole data. 

We analyze with MCMC and compare three classes of cosmological models: 

• The usual slow-roll ACDM, the ACDM+r and the ACDM+r on Cbni models. BNI stands for Binomial New 
Inflation. In this last model we enforce the theoretical functional relation (denoted Cbni) between n s and r 
valid in BNI. (We call ACDM+r on Cbni the usual ACDM+r model constrained on the curve Cbni)- 

• The slow- roll ACDM on Cbni model with a sharp cut for k < k\. 

• The ACDM on Cbni model including both fast and slow-roll stages. 

The MCMC analysis of the WMAP and SDSS data favours a double- well, spontaneously broken symmetric potential 
for the inflaton in new inflation [15J 



Here N ~ 60 is the number of efolds since the cosmologically relevant modes exit the horizon till the end of inflation 
and Mpi = 2.4 x 10 18 GeV is the Planck mass. MCMC yields for the dimensionless quartic coupling y ~ 1.32 and 



for the inflation energy scale M and the inflaton mass scale m, respectively. 

We modified the CosmoMC code introducing the fast roll transfer function D-ji(k) in the primordial power spectrum 
according to eq. (|3.3ip . 

We repeated the same analysis with the WMAP-5, ACBAR-2007 and SDSS data, this time setting N = 60, with 
no statistically significant change. 

Our fits imposing Cbni predict non-zero lower bounds on r: at 95% CL, we find that r > 0.023 when no cutoff is 
introduced, while r > 0.018 when either the sharp cutoff or the fastroll D(k) are introduced. The best fit values of 
the other cosmological parameters remain practically unchanged as compared to ACDM. Similarly their marginalized 
probability distributions are almost unchanged, with the natural exception of n s , which in BNI has a theoretical upper 
limit [see eq. (|5 ■ - 

We observe that the oscillatory form of the fastroll transfer function D-jz(k), by depressing as well as enhancing 
the primordial power spectrum at long wavelengths, leads also to new superimposed oscillatory corrections on 
the low multipoles. As far as fitting to current data is concerned, such corrections are more effective than the pure 
reduction caused by a sharp cutoff. The fast-roll oscillations yield better gains in likelihood than the sharpcut case. 

We display the best fit for the Cj T , Cj E and Cf E multipoles compared to the experimental data at low I. One 
can observe that for £ — 2 and £ — 3 fast-roll and sharpcut models yield rather similar results (and better than the 
ACDM+r model) while for I = 4 fast-roll produces for Cj E a value closer to WMAP-3 data than sharpcut. For C EE 
both fastroll and sharpcut models produce a depression of the low multipoles including the EE quadrupole. 

We summarize in the Appendix the numerical code used by us in the simulations. 

We display the real space two point TT-correlator C TT {9) for purely slow-roll ACDM, sharpcut and fast-roll ACDM 
models. The purely slow roll ACDM correlator differs from the two others only for large angles 9 > 1. Since all /-modes 
besides the lowest ones are practically identical in the three cases, this shows how important are the low multipoles 
in the large angle correlations. 

We get the following picture of the inflationary universe explaining the quadrupole suppression from the effective 
(Ginsburg-Landau) theory of inflation combined with MCMC simulations of CMB+LSS data. A fast-roll stage lasting 
about one efold is followed by a slow-roll stage lasting ~ 65 efolds. We have the radiation dominated era after these 
~ 65 + 1 = 66 efolds of inflation. The quadrupole modes exit the horizon during the fast-roll stage about 0.4 of an 




The quartic coupling in the effective theory of inflation [5J is given by 




M = 0.57 x 10 lb GeV 



m = 1.34 x 10" GeV 
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cfold after the beginning of inflation and is therefore suppressed compared with the modes exiting later the horizon 
during the slow-roll stage. 

The fast-roll stage explains the quadrupole suppression and fixes the total number of efolds of inflation. The 
fact that the quadrupole mode £;q exits the horizon before the slow-roll stage implies an upper bound in the total 
number of efolds N to t during inflation. Combining this with estimates during the radiation dominated era we obtain 
TV > 56, Ntot ~ 66, the upper bound Ntot < 82 and the lower bound Not > 62. 

Our MCMC simulations give good fits for TV = 50 and TV = 60. The bound TV > 56 therefore favours TV ~ 60 which 
implies TV tot ~ 66 and H ~ 3 x 10 10 GeV by the end of inflation. 

Changing TV from 50 to 60 does not affect significatively the MCMC fits we present in this paper. This is partially 
due to the fact that a change on y can partially compensate a change on TV. More importantly, a 20% change in TV 
may affect the fit of k\ by a similar amount, leaving unchanged its scale, which is of the order of the inverse Hubble 
scale. 

Another hint to increase TV above 50 comes from WMAP-5 that gives a larger n s value and using the theoretical 
upper limit for n s [3, EH]: n s < 1 — 1 ' 9 ^ 6 '" , which gives n s < 0.9679 ... for TV = 60. This value is compatible with 
the n s value from WMAP5+BAO+SN and no running [!|. 



II. THE EFFECTIVE THEORY OF INFLATION. FAST AND SLOW ROLL REGIMES. 

The inflaton potential V(cp) must be a slowly varying function of ip in order to permit a slow-roll solution for the 
inflaton field which guarantees a large enough total number of efolds > 62. Such value is necessary to solve the 
horizon, flatness and entropy problems. 

As discussed in ref. Q , the inflaton potential should have the universal form 

V(<p)=NM i w( X ), (2.1) 

where x is a dimensionless, slowly varying field 

x = 7W^r (2 - 2) 

and M is the energy scale of inflation, TV ~ 60 is the number of efolds since the cosmologically relevant modes exit 
the horizon till the end of inflation. 

The energy scale M of inflation is determined by the amplitude of the observed CMB anisotropy, which implies 
M ~ 0.7 x 10 16 GeV. That is, M <C Mpi, which ensures the consistency of the effective theory of inflation. 

The dynamics of the rescaled field x exhibits the slow time evolution in terms of the stretched dimensionless time 
variable, 

tM2 Ws (2.3) 



Mpi VTV VN M 2 

The rescaled variables x an d t change slowly with time. A large change in the field amplitude ip results in a small 
change in the x amplitude, a change in <p ~ Mpi results in a x change ~ 1/VTV- The form of the potential, eq. (l2.1|l . 
the rescaled dimensionless inflaton field eq. ()2.2j) and the time variable r make manifest the slow- roll expansion as a 
consistent systematic expansion in powers of 1/TV [5J]. 

We can choose |u>"(0)| = 1 without loosing generality. Then, the inflaton mass scale at zero field is given by a 
see-saw formula 

M 4 M 2 

m 2 = \V"(<p = 0)| = — 2- , M~0.7x 10 16 GeV , m=— 2.0 x 10 13 GeV . (2.4) 

Mpi Mpi 

The Hubble parameter when the cosmologically relevant modes exit the horizon is given by 

H = VN mH - 5 to - 1.0 x 10 14 GeV , (2.5) 

where H ~ 1. As a result, m <C M and H <C Mp\. 
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The energy density and the pressure in terms of the dimensionless rescaled held x an d the slow time variable r take 
the form, 

The equations of motion in the same dimensionless variables become 

1 



27f(*) +l ° (x) 



The slow-roll approximation follows by neglecting the terms in eas. (|2.7p . Both w(x) and 7i(r) are of order iV° for 
large N. Both equations make manifest the slow roll expansion as an expansion in 1/N. 

The number of e-folds N[x] since the field \ exits the horizon till the end of inflation (where x takes the value Xend) 
can be computed in close form from eqs. (|2.7p in the slow-roll approximation (that is, neglecting 1/N corrections): 

_ /*- -(x) dx 1 (2 . 8) 



7 X «/(x) 

where we choose N = -ZV[x]- Therefore, eq. (|2.8[) determines x &t horizon exit as a function of the couplings in the 
inflaton potential w(x) : 

^Ld X =l. (2.9) 
w'{x) 

Inflation ends after a finite number of efolds provided 

w(Xend) = w'{Xend) = . (2-10) 

So, this condition is enforced in the inflationary potentials. 

There are two generic inflationary regimes: slow-roll and fast-roll depending on whether [7( 

1 fdx\ 2 

-^—^ ( — ) <C w(x) ■ slow — roll regime 

w(x) ■ fast — roll regime . (2-11) 



2 N V dr 



1 fd X " 2 



2 N \dr 

Both regimes appear in all inflationary models in the class ea. (|2.1|l . Fast-roll clearly corresponds to generic initial 
conditions for the inflaton field. The fast-roll stage turns to be very short and is generically followed by the slow-roll 
stage 0- 

For the quartic degree potentials V(tp), the main two families are: 

(a) discrete symmetry (ip — > —tp) breaking potentials (so-called new, or small-field, inflation) 

. , A / 9 m 2 \ 2 m 2 9 A , m 4 . „. 

^ = 4 V Xj + 4A J (2 - 12) 

(b) unbroken symmetry potentials (chaotic, or large-field, inflation), 

%)=4^ 2+ i»' 4 - (2 - i3) 

The corresponding dimensionless potentials w(x) take the form 

w(x) = ^(x 2 --) =-\x 2 + ^X i + ~ for new inflation (2.14) 
32 V y ) 2 32 y 
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and 



1 V 

w(x) — 2 X 2 + 22 ^ or cna °ti c inflation 



(2.15) 



where the coupling j/ is of order one and 



A = 



y ( M 



N \Mpi 



1CT 



In new inflation the inflaton starts near the local maximum \ = and keeps rolling down the potential hill till it 
reaches the absolute minimum \ = yf ■ ^ n general, the initial kinetic energy may be of the same order of magnitude 
as the initial potential energy of the inflaton which defines fast-roll inflation. That is, in general the initial states are 
not slow-roll. 

By numerically solving eqs. (|2.7p we find that the fast-roll initial stage of the inflaton becomes very soon a slow-roll 
stage 0- This is a general property and implies that the slow-roll regime is an atractor for this dynamical system 
[|[. We see a de Sitter-like expansion during the slow-roll stage r < 3 during which the Hubble parameter decreases 
slowly and monotonically. 

An initial state for the inflaton (inflaton classical dynamics) with approximate equipartition between kinetic and 
potential energies is a more general initialization of cosmological dynamics in the effective field theory than slow 
roll which requires that the inflaton kinetic energy is much smaller than its potential energy. The most generic 
initialization of the inflaton dynamics in the effective field theory leads to a fast roll stage followed by slow roll 
inflation 0]. 

The total number of efolds of inflation is determined by the initial conditions for the inflaton field: x(0) ~ x(0) = 
0(1). Varying these initial conditions the total number of efolds of inflation sweeps a wide range of efold values 
showing the flexibility of the inflationary model. 

We have carried out analogous numerical studies in scenarios of chaotic inflation with similar results: if the initial 
kinetic energy of the inflaton is of the same order as the potential energy a fast roll stage is always present. The 
evolution of the potentials Viz(r]) and Vt(v) felt by the curvature and tensor perturbations are similar to those for 
new inflation and they are always attractive during the fast roll stage (see below). 




FIG. 1: We plot here e v vs. In a during the fast roll stage and the beginning of slow-roll for new inflation with y — 1.322. We 
define as the end of fast-roll the point where e v — -h = 0.02. This gives here lna^ = 1.091. Namely, fast-roll ends one efold 
after the beginning of inflation. 



III. THE EFFECT OF FAST-ROLL ON THE INFLATIONARY FLUCTUATIONS. 



The inflationary scenario features scalar curvature fluctuations determined by a gauge invariant combination of the 
inflaton field and metric fluctuations. They also feature tensor fluctuations (gravitational waves). 
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It is convenient to introduce the gauge invariant potential Q , 

u(x, t) = — z lZ(x, t) . (3-1) 
where lZ(x,t) stands for the gauge invariant curvature perturbation of the comoving hypersurfaces and 



z = a(t) 



H 



(3.2) 



The gauge invariant curvature field u(x,t) expanded in terms of conformal time mode functions and creation and 
annihilation operators take the form [2| 



U(X,T)) 



d 3 k 
(2 7T)i 



Q»(fe) S n (k; rj) e lk x + o^(k) S^(k; V ) e 



ik-x 



(3.3) 



where the operators obey canonical commutation relations 

a n (k), <4(fc')l =5 {3 Hk-k') . 
The vacuum state is annihilated by the operators ttn(k) and the mode functions obey the equations of motion 0], 



d 2 k2 1 d 2 z 
drf z drf 



S n (k;r,) = 0. 



Here rj stands for the conformal time 



V 



dt 



Eq. (|3.4I) is a Schrodinger-type differential equation in the variable r\. The potential felt by the fluctuations 

1 d 2 z 
z dr\ 2 



W M = z — 



(3.4) 



(3.5) 



(3.6) 



can be expressed in terms of the inflaton potential and its derivatives. From eqs. (|3.2[) and (|3.6[) and using the inflation 
equations of motion p.7p . the potential Wn{rf) can be written as @ 



2-7e v + 2e 2 v 



e v V 



M Pl H 2 

where we take for the sign of the square root J(Z the sign of ip and 



Vv{3 - e v ) 



1 V 2 M 2 V"{ V ) 



2M 2 pl H 2 



(3.7) 



(3.8) 



e v and r) v are the known slow-roll parameters 0. Notice that eqs. (|3.7[) - p.8[) are exact (no slow-roll approximation). 
In terms of the dimensionless variables eqs. (|2.1|) - (|2.3|) we obtain for the potential Wm{rj), 



Wn{ri) =a J (r)) Ti m 2 N 



while the parameters e v and r\ v take the form 



2 - 7 e v + 2 e 2 v 



n n 2 



r) v (3 - e v ) 



1 1 (i X 



" 2 N H 2 \dr 
In the slow-roll regime they can be approximated as 

2 



Vv = TF 



1 w"( X ) 
N w( X ) 



1 



w'(x) 



N 2 



Ol i ] = O I - 



N 



_ 1 w"{ X ) 
VV ~ N w( X ) 



(3.9) 



(3.10) 



(3.11) 
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We explicitly see that the parameters e v and r\ v are suppressed by powers of l/N in the slow-roll regime. This result 
is valid for all models in the class defined by eq. (|2.ip regardless of the precise form of w(x)- 

Tensor perturbations (gravitational waves) are gauge invariant. The corresponding quantum fields (gravitons) are 
written as 



hUx,rj) = 



a(rj) Mpt 



E 

A=x,+ 



d 3 k 

(2^1 



4(A,fe) 



flA.fe S T (k,r)) + e lkx a\ k S^(k,rj) 



(3.12) 



where A labels the two standard transverse and traceless polarizations x and +. The operators a\.k, & obey 
canonical commutation relations, and £ij(A, k) are the two independent symmetric and traceless-transverse tensors 
constructed from the two independent polarization vectors transverse to k, chosen to be real and normalized such 
that e}(A,fc) e{{X',k) = 8\ 6 X ,\>. 

The mode functions St (&;??) obey the differential equation 0, [HQ 



S T (k;rj) 



a"(v) 



a(rj) 



S T {k;r 1 ) = 



(3.13) 



That is, for both scalar curvature and tensor equations we have the equation 

r d 2 



dvj 2 



W(rf) 



S(k;r,) = 



(3.14) 



where for scalar curvature perturbations Wn(rj) is given by eq. (|3.6|) and for tensor perturbations Wt(t)) is 

a"(v) 



W T (V) 



a{rj) 



It is convenient to explicitly separate the behavior of W(ji) during the slow roll stage by writing 

,,2 1 



(3.15) 



(3.16) 



W(jj) = V(77) + , 

where the potential V{rj) is the fast-roll part and, 

fv-ji = | + 3 e v — rj v + O (wz) for curvature perturbations 
vt = | + e v + O (j^r) for tensor perturbations . 

e v and T) v are given by eqs. (|3.8p . (|3.10p . 

The potential V(t]) is localized in the fast roll stage prior to slow roll (during which cosmologically relevant modes 
cross out of the Hubble radius), V(j]) vanishes during slow-roll. In terms of the potential V{rf) the equations for the 
quantum fluctuations read, 



drf 



T 



V(r?) 



S(k; V ) = 



(3.17) 



A. The Fluctuations during the slow-roll stage 

The slow roll dynamics acts through the term 

v 2 - 1/4 
rf 2 



which is a repulsive centrifugal barrier. 
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During the slow roll stage V(ry) is negligeable and the mode equations simplify to 



drj 2 



,2 i 



S sr (k, rj) = 



(3.18) 



To leading order in slow roll, v is constant and for general initial conditions the solution of eq. (|3.18p is, 

S sr (k; yy) = A(k) g u (k; rj) + B(k) [g v {k; rj)]* , (3.19) 

where 



are Hankel functions. These solutions are normalized so that their Wronskian is given by 
W[g v (k; rj), #*(fc; rj)] = g' v (k; rj) g* v {k; rj) - g v (k; rj) [g' v {k; rj)]* = -i . 



(3.20) 



(3.21) 



The mode functions and coefficients A(k), B{k) will feature a subscript index 1Z, T, for curvature or tensor pertur- 
bations, respectively. 

For wavevectors deep inside the Hubble radius 3> 1, the mode functions have the asymptotic behavior 



9v{k;rj) = - 
while for 77 — ^ — the mode functions behave as: 

?7^o" r(i/) 



1 



-ikr/ 



<2k 



V2tt k \ikij 



(3.22) 



(3.23) 



In particular, in the scale invariant case v = I which is the leading order in the slow roll expansion, the mode functions 
eqs. (|3.20j) simplify to 



9i(k;v) 



-ikr/ 



2k 



krj 



(3.24) 



B. The Fluctuations during the earlier fast-roll stage 



The mode equation (|3.17[) can be written as an integral equation. We choose as initial condition the usual Bunch- 
Davies asymptotic condition 



S(k; rj — > —00) = g v {k; rj — > —00) 



-ikr) 



2k 



(3.25) 



We formally consider here inflation and the conformal time starting at 77 = —00. However, it is natural to consider 
that the inflationary evolution of the universe starts at some negative value rji < fj, where 77 is the conformal time 
when fast roll ends and slow roll begins. 

The mode equation (|3.17[) can be written as an integral equation including the Bunch-Davies initial condition 

S(k;ri)=g„(k; V )+ig u (k;ri) f g*.(k; rf) V(r/) S(k; tf) drf - i g* v {k; rj) f 9v{k; n') V(rf) S(k; r,') drf . (3.26) 



where for simplicity we set 77, = —00. 

Since V{rj) vanishes for 77 > 77, the mode functions S(k; rj) for rj > fj can be written as linear combinations of the 
mode functions g u (k] rj) and g*(k; 77), 



S{k;rf) = A{k) g v {k;rj) + B{k) g* v {k;rj) , 77 > 77 



(3.27) 
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where the coefficients A(k) and B{k) can be read from ea. (|3.26p . 



A(k) 



1+i 



9t(k;v) V(v) S(k;rj) dr/ 



B(k) = -i f g v (k;v)V(ri) S(k; V ) d v 

J — oo 



(3.28) 



(3.29) 



The coefficients A(k) and B(k) are therefore calculated from the dynamics before slow roll, that is, during fast-roll, 
[recall that V {rj) = for rj > fj during slow roll.] 

The constancy of the Wronskian W[S(k; rj), S*(k; rj)] = —i and eqs. (|3.21|) . (|3.27p imply the constraint, 

\A(k)\ 2 -\B(k)\ 2 = l . 
This relation permits to represent the coefficients A(k); B(k) as |6j 

A(k) = ^1 + N{k) e l ° A(k) ■ B{k) = ^N{k) e l0B(k) , (3.30) 

where N(k), &A,B(k) are real. 

Starting with Bunch-Davies initial conditions for rj — > — oo, the action of the fast-roll potential V(j]) generates 
a mixture (Bogoliubov transformation) of the two linearly independent mode functions <?i/(fc;f?) and g*{k;rj), which 
result in the mode functions S{k;rj) eq. (l3.27p for rj > fj when the fast roll potential V(r?) vanishes. This is clearly 
equivalent to starting the evolution of the fluctuations at the beginning of slow roll r\ = fj with initial conditions 
defined by the Bogoliubov coefficients A(k) and B(k) given by eq. (|3.29[) as stressed in ref. Q- 

As shown in ref. [7j the power spectrum of curvature and tensor perturbations for the general fluctuations eq. (|3.27|) 
takes the form, 



Pn(k) 



,o- k 3 



2 7T 2 



P T (k) " = L 



o- k 3 



2 n 2 



z{rf) 
S T (k;n) 



P%(k) 1 + D n (k) 



= Pf(k) l + D T (k) 



(3.31) 



Here Dn(k) and Dr{k) are the transfer functions for the initial conditions of curvature and tensor perturbations 
introduced in ref. [6|: 



D(k) = 2 \B(k)\ 2 - 2 Re [A(k) B*(k) 



2v-31 



2 N(k) - 2 y/N(k)[l+N(k)] 



(3.32) 



where one uses either r or t quantities and 9k = ds{k) — $A{k). 
The standard slow roll power spectra are given by @ : 



Pr(k) 



k 

2h 



H 2 



r 3 2 e v M 2 



A 2 



n 



Al 



ko 



pi 
2e,. 



A 2 



A 2 



16 £ v 



(3.33) 



These spectra are modified by the fast-roll stage as displayed in eq. (|3.31[) . The scale ko is a reference or pivot scale, 
for example WMAP takes k = 0.002 Mpc" 1 and CosmoMC, k = 0.050 Mpc" 1 (see section ITVT) . 

The integral equation (13.26|) can be solved iteratively in a perturbative expansion if the potential V{nj) is small 
when compared to 

,2 - 2 -V4 



which is indeed true in this case. Then, we can use for the coefficients A[k), B(k) the first approximation obtained 
by replacing S(k; rj') by g v (k; rj') in the integrals eqs. (|3.28|) - (|3.29p . This is the Born approximation, in which 



A(k) 



1 + i 



V(ri)\g„(k;r,)\ 2 dr, , B(k) 



V{v)gl(k;v) drj . 



(3.34) 
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The transfer function of initial conditions given by eq. (|3.32[) can be computed in the Born approximation, which is 
indeed appropriate in this situation. By using eqs. (|3.34[) for the Bogoliubov coefficients A(k) and B(k) to dominant 
order in 1/N, that is v = 3/2 [eq. (|3. 16|> ] . D(k) is given by, 



D(k) 



dr, V(t?) 



sin(2 k ri) 1 - 



The potential V{rj) is obtained from eq. ()3.15p as 



V(rj) = W(r)) 



k 2 rj 2 



v 2 - 1/4 



k r\ 



cos(2 k rj) 



(3.35) 



To explicitly compute V-jiivj) as a function of rj for the curvature fluctuations we solve numerically the equations of 
motion (|2.7p for new inflation [eq. (|2.14[l ] and insert the solution for the inflaton xiv) m eqs. (|3.9[) - (|3.10[) . No large N 
approximation is used in this numerical calculation since we cover in the evolution the fast-roll region where slow-roll 
obviously does not apply. 

We plot in fig. [2] V-jziv) vs - V f° r new inflation [eq. (|2.14[) ] for the coupling y = 1.322 and a total number of efolds 
equal to sixty. We choose here the initial values of \ an d X such that their initial kinetic and potential energies are 
equal. We see that the potential V-jiixi) is attractive in the fast-roll stage and asymptotically vanishes by the end of 
fast roll i] — -0.04. 



-50 



-100 



-150 



-200 



-250 




-0.25 



FIG. 2: The potential Viz(rj) vs. rj for new inflation with y — 1.322. Viz{v) i s attractive during fast-roll and vanishes by the 
end of fast roll (rj ~ -0.04). 



We obtain the transfer function D-jz(k) by inserting V-ji{rj) m f° eq. (|3.35|) and computing the integral over rj numeri- 
cally. In fig. [3] we plot D-ji(k) vs. k/m for new inflation [eq. (|2.14p ] and ten different couplings 0.00536 < y < 1.498 with 
a total number of efolds equal to sixty. We see that Dn(k) oscillates around zero and therefore produces suppres- 
sions as well as enhancements in the primordial power spectrum [see eq. (|3.31|) ]. D-n(k) vanishes asymptotically 
for large k as expected. 

The first peak in Dn{k) is clearly its dominant feature. The k of this peak corresponds to k- modes which are today 
horizon size and affect the lowest CMB multipoles (see below and table 2) 0,0]. 

For small k the Born approximation to Dn{k) yields large negative values indicating that this approximation cannot 
be used in this particular small k regime. We introduce the scale k\ by the condition Dn{kx) = — 1 and then just 
take Dn{k) — — 1 for k < k\. This corresponds to vanishing primordial power for the lowest values of k [see fig. [3]. 

^From fig. [3] we also see that the plots of D-ji(k) for different couplings follow from each other almost entirely 
by changing the scale in the variable k as summarized by eq. (|3.36p . Indeed, the characteristic scale k\ plays here a 
further important role. 

Analysing V-jiiv) and Dn(k) for different couplings y we find that they scale with k%. Namely, 



(3.36) 
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FIG. 3: D-Ti(k) vs. k/m for new inflation and ten different couplings 0.00536 < y < 1.498. We see that the plots of D-n(k) for 
different couplings follow from each other by changing the scale in the variable k as summarized by eq. (|3.36[ l. 
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FIG. 4: fci vs. y for new inflation. 



where Q(x) and ^(x) are universal functions. That is, Q{x) and ^{x) do not depend on the coupling y while 
ki = ki(y). We display k\ vs. y in fig. Q] 

We obtain the function Q{x) from ea. (|3.36|) as, 



(3.37) 



We plot Q(x) in fig. [5] as follows from the r. h. s. of eq. (|3.37|) for ten different values of y. We see that all the curves 
collapse on a common curve proving the validity of the quasi-scaling properties eq. (|3.36[) . 
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FIG. 5: Q(x) for the ten values of y of fig. [3] according to eq. (|3.37| l. All curves collapse to a common one proving the scaling 
properties eq. (|3.36[) . 



IV. MCMC ANALYSIS OF CMB AND LSS DATA INCLUDING THE EARLY FAST-ROLL 

INFLATIONARY STAGE 

In order to test the theoretical quadrupole depression predicted by fast-roll inflationary stage against the current 
experimental data we performed a Monte Carlo Markov Chains (MCMC) analysis of the commonly available CMB 
and LSS data using the CosmoMC program 

For LSS we considered SDSS (DR4). For CMB we first considered the three-years WMAP data (with the second 
release of WMAP likelihood code) and small scale data (ACBAR-2003, CBI2, BOOMERANG03). While this work 
was in progress the five-years WMAP data were released, and we repeated our MCMC analysis almost completely 
with these new data, using also the newer 2007 ACBAR release. Actually WMAP-3 or WMAP-5 provide by far the 
dominant contribution and small scale experiments have very little relevance for the quadruple depression issue. 

In all our MCMC runs we have not marginalized over the SZ (Sunayev-Zel'dovich) amplitude and have not included 
non-linear effects in the evolution of the matter spectrum. The relative corrections are in any case not significant 
Hi EH i especially in the present context. 

CosmoMC is a publicly available open-source FORTRAN package that performs MCMC analysis over the parameter 
space of the Standard Cosmological model and variations thereof. The main observables in this approach are the cor- 
relations among the CMB anisotropics and in particular: the TT (Temperature- Temperature) , the TE (Temperature- 
Erodes), the EE (E_modes-E_modes) and the BB (B_modes-B .modes) correlation multipoles (E_modes and B_modes 
are special modes of the CMB polarization) . These multipoles can be numerically calculated with very good accuracy, 
as functions of the cosmological parameters, from the primordial power spectrum through programs such as CAMB 
(included in CosmoMC). On the other side, experimental data provide a likelihood distribution for multipoles, which 
is then turned into a likelihood for the cosmological parameters through the MCMC method. We modified the 
CosmoMC code introducing the transfer function D-r (k) in the primordial power spectrum according to eq. (|3.3ip . 

We ran CosmoMC on pc clusters with Message Passing Interface (MPI), producing from 10 to 24 parallel chains, 
with the 'R-l' stopping criterion (which looks at the fluctuations among parallel chains to decide when to stop the 
run) set equal to 0.03. The statistical converge was also verified a posteriori with the help of the getdist program of 
CosmoMC. 

The preferred reference model for slow-roll inflation cosmology is the ACDM+r model, that is the standard ACDM 
model, which has six parameters 1 , augmented by the tensor-scalar ratio r. Indeed, the current experimental accuracy 
provides sensible bounds only for the first order parameters e v and rj v , through their standard relation to the scalar 
spectral index n s and the ratio r : n s — 1 = 2 r\ v — 6 e v , r — 16 e v . Specific slow-roll scenarios, such as those based on 



1 we use the standard ones of CosmoMC, that is the baryonic matter fraction Wb, the dark matter fraction uj c , the optical depth r, 
the ratio of the (approximate) sound horizon to the angular diameter distance 8, the primordial supcrhorizon power in the curvature 
perturbation at 0.05 Mpc - x , A s and the corresponding scalar spectral index n s 
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new (small-field) or chaotic (large-field) inflation, connect in a model-dependent way e v and r) v to free parameters in 
the inflaton potential and thus typically lead to specific theoretical constraints in the (n s , r) plane 15]. 

We point out that we used the default CosmoMC pivot scale ko = 0.05 Mpc -1 rather than the customary WMAP 
choice of fco = 0.002 Mpc" 1 . As evident from eq. (|3 . 33[) this leads to a small difference with respect to the WMAP 
choice in the definition itself of the tensor-scalar ratio r. In particular, the CosmoMC r is roughly 10% larger than 
the WMAP one. 



A. MCMC analysis without quadrupole suppression: Dn(k) = 0. 

Let us present our MCMC analysis with the standard slow-roll primordial power eq. (l3.33j) . That is, without 
including the early fast roll stage and therefore vanishing transfer function Dfi(k). 

For instance, in the simplest binomial realization of new inflation described by the inflation potential of eq. (|2.12|) 
or eq. (|2.14p . n s and r are constrained to the curve Cbni (BNI stands for Binomial New Inflation) parametrized by 
the quartic coupling y as (lEj : 

ns = l 'NjY^zf' r= ~(T^F' 2/ = ^-l-log^ , ^gX 2 , 0<z<l. (4.1) 

This situation is clearly displayed in fig. [6] in which the curve Cbni, for the two choices TV = 50 and N = 60, is drawn 
over the contour plot of the likelihood distribution for n s and r in the ACDM+r model obtained with CosmoMC, 
using the WMAP-3, small-scale CMB and SDSS data. Practically, the same contour plot applies when WMAP-5 and 
ACBAR-2007 are used. 

The likelihood L, as function of the whole set of parameters, provides a quantitative measure of the power of a given 
model to fit the multipoles CJ. As customary, we set — 21ogL = Xl> although it is well known that, due particularly 
to cosmic variance, the shape of L, as function of the C7, is not Gaussian especially for low I. 

Now, as evident from eq. (|4.ip and fig. [51 one could expect from the ACDM model constrained to Cbni a fit to the 
data not as good as in the ACDM+r model since the current data seem to favor smaller values for r. Indeed we find 

minx! (ACDM + r on Cbni) - minx! (ACDM + r ) - 0- 4 • (4-2) 

This result was obtained for N = 50 by direct minimization of x! hi the neighbourhood of Cbni, using the data of a 
large collection of long chain runs (with a grandtotal of almost two million steps) for the ACDM+r model with the 
WMAP-3, small-scale CMB and LSS data. The flat priors on the cosmological parameters were the standard ones of 
CosmoMC, that is 

0.005 < u b < 0.1 , 0.01 < uj c < 0.99 , 0.5 < 6 < 10 
0.01 < t < 0.8 , 2.7 < log(10 10 yl s ) < 4 , 0.5 < n s < 1.5 

while for the tensor-scalar ratio we imposed as prior 

< r < 0.35 . 

We repeated the same analysis with the WMAP-5, ACBAR-2007 and SDSS data, this time setting N = 60, with no 
statistically significant change. 

Another approach, that unlike the direct minimization of x! over Cbni does take advantage of the explicit analytic 
parametrizations in eq. (|4.ip , is to use the single variable z as MCMC parameter, instead of the constrained pair 
(n s r), with a flat prior over all the allowed range < z < 1. Let us call ACDMz|Cbni the 6-parameter model ACDM 
constrained on Cbni using the variable z. Then, we find that taking into account the natural fluctuations due to 
the large number of data (which make the likelihood landscape over the MCMC parametrs quite complex) and the 
various approximations and numerical errors in the theoretical calculation of the multipoles, the increasing in x! due 
to the Cbni eq. (|4.ip constraint compared to the ACDM model essentially vanishes [see Table I below] . 

For completeness and reference, we report in Table I our best fit (or most likely) values for the MCMC cosmological 
parameters, as well as the absolute value of our best likelihoods, which of course depend on the specific datasets used, 
that is WMAP-3, small scale CMB and SDSS. We report in the first line of Table I our best fit for the standard 
ACDM model, which has six free parameters since r is set to zero by fiat. 

It should be noted that the likelihoods difference between ACDM+r and its direct restriction to Cbni, (that is 
ACDM +r on Cbni), is mostly due to the SDSS data, which together with the CMB data place an upper bound 
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FIG. 6: Binomial New Inflation compared to ACDM+r model in the (n s , r) plane. The color-filled areas correspond to 
12%, 27%, 45%, 68% and 95% confidence levels for ACDM+r according to WMAP-3, small scale CMB and SDSS data. 
Cbni is the solid red curve for N = 50 or the dashed magenta curve for N — 60. The white dots corresponds to the values 
0.01 + 0.11 * n, n = 0, 1, . . . , 9, of the variable z in eq. (|4.ip . starting from the leftmost ones. The quartic coupling y instead 
increases monotonically starting from the uppermost dots, which corresponds to the free-field, purely quadratic inflaton potential 
y — 0. We see that very small values of r are excluded since they correspond to n s < 0.92 outside the 95% confidence level 
contour. 
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10 9 A S 


n s 


r 


xi/2 


ACDM 


0.224 


0.106 


1.041 


0.886 


2.072 


0.959 





2713.906 


ACDM+r 


0.224 


0.107 


1.042 


0.831 


2.054 


0.960 


0.009 


2713.972 


ACDM+r on Cbni 


0.223 


0.106 


1.040 


0.848 


2.047 


0.956 


0.059 


2714.166 


ACDMz|C B ni 


0.222 


0.107 


1.041 


0.877 


2.065 


0.958 


0.069 


2713.918 



TABLE I: Best fit values for the MCMC cosmological parameters without quadrupole suppression, using WMAP-3, small-scale 
CMB and SDDS. Cbni means the curve on which n s and r are constrained in Binomial New Inflation (BNI), eq. (|4.1|) with 
N = 50. ACDM+r on Cbni means the ACDM+r model constrained on Cbni- ACDMz|Cbni denotes the ACDM model 
constrained on Cbni using the single variable z eq. (|4.1[) as MCMC variable instead of the constrained pair (n s ,r). 

on r twice more stringent than WMAP-3 alone. Indeed, when only the WMAP-3 data were used, we verified that 
no significant likelihoods difference was exhibited. This situation changes slightly when WMAP-5 is used, with xi 
increasing approximately by 0.2 from ACDM+r to ACDM+r on Cbni, since WMAP-5 alone puts a tighter bound on 
r than WMAP-3 alone (0.43 vs. 0.65 at 95% CL [§]). 

It is evident that, as far as most likely values of the cosmological parameters are concerned, the fit with the constraint 
Cbni included, either with or without z as MCMC parameter, does not determine any statistically significant change, 
except of course for n s and r themselves. In particular, with respect to the ACDM + r results, the most likely value 
of n s is practically unchanged, while that of r changes from values of order 10~ 2 (or just in ACDM) to values such 
as 0.059 and 0.069 (see Table I). 

Concerning marginalized distributions, we find no significant changes but for n s and r. These results are very 
close to those in ref. [l5l |. where trinomial new inflation, with a possibly asymmetric potential, was considered. In 
particular, the marginalized distribution for r shows a broad but clear peak centered near the most likely value as in 
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[l5l |. In the present context of binomial new inflation, we find that r = 0.0891q'q5 4 with r > 0.023 at 95% CL. 

All together, these results show that, as far as pure data fitting is concerned, the 6-parameter ACDMz|Cbni model 
is just as good as the standard 6-parameter ACDM model. More generally speaking, we may say that current CMB 
and LSS data together, without any theoretical constraint, put only an upper bound on r (namely r < 0.20 with 95% 
confidence level in the most recent WMAP-5 analysis [9j). Therefore, any inflation-based 6-parameter model (as the 
ACDMz|Cbni model, for instance) predicting a value of r well below 0.2 is as likely as the ACDM model itself. This 
means that the theoretical grounds of a given model take a more important role in the analysis and interpretation 
of the CMB and LSS data. For instance, from an inflationary viewpoint, the choice that r vanishes exactly appears 
unlikely and unphysical. Notice that n s — 1 = r = corresponds to a singular and critical (massless) limit where 
the inflaton potential vanishes [lj|, while The MCMC analysis for both the binomial and trinomial new inflationary 
models yield lower bounds for r. 

In order to asses the statistical relevance of the quadrupole suppression, we determine in the best fit ACDM model, 
the probability that there is at least one multipole, regardless of £, smaller that 20% of the theoretical mean value. We 
obtained 0.06126 for such probability. Thus, in the ACDM, the observed quadrupole realizes a rather unlikely event 
which has only a 6% probability. Therefore, it makes sense to search for a cosmological explanation of the quadrupole 
depression beyond the ACDM model. 



B. MCMC analysis including the quadrupole suppression: Djz(k) ^ 0. 

Let us now further develop this argument by considering the quadruple depression, avoiding the a priori dismissal 
based on the simple invocation of cosmic variance or experimental inaccuracy. In the standard ACDM model the 
simplest, purely phenomenological way to decrease the low multipoles is to introduce a infrared sharp cut in the 
primordial power spectrum of the curvature fluctuations. That is, one assumes that Pn(k) = for k < ki and treats 
ki as a new MCMC parameter to be fitted against the data. It is actually not necessary to include also a cut on 
the tensor power spectrum, since it would lead to changes certainly not appreciable within the current experimental 
accuracy. 

With this procedure we obtained, using either the WMAP-3 data alone or both CMB and LSS data: 

minx! (ACDM + sharpcut) - minx! (ACDM) - -1-4 

This result is slightly better than the one reported in ref. 8|, but still the likelihood gain hardly compensates the 
price of a new parameter, especially because its nature appears quite ad hoc. In fig. [7] we plot the marginalized 
probabilities and mean likelihoods of the seven MCMC parameters plus other standard derived parameters in the 
CMB+LSS case. In the WMAP-3 alone case these plots are almost identical. There are no significant changes from 
ACDM to ACDM+sharpcut in their common parameters, in either most likely values or marginalized distributions. 
The distribution of the new cutoff parameter k\ shows a well defined peak centered on its most likely value (ML), 
which corresponds to today's physical wavelength 



0.291 (Gpc)- 1 (WMAP-3 only) 
0.272 (Gpc)- 1 (CMB+LSS) 



( fc i)ML = \ :::: tz' ; z ( acdm + ^cut) , 



that is of the order of today's inverse Hubble radius, as expected. 

Introducing the infrared sharp cut on Pji(k) in the ACDMz|Cbni model we find sizably different gains 



minx!(ACDMz|C B Ni + sharpcut) - minx!(ACDMz|C] 



BNI 




(WMAP-3 only) 
(CMB+LSS) 



As before, the difference is due to the tighter bound on r due to the inclusion of the SDSS data. In fact, the most likely 
values (ML) of k\ and r corresponding to min x!(ACDMz|Cbni + sharpcut) are given in Table II. The marginalized 
probabilities in the (r, fci) plane (converting to a flat prior on r) are shown in the two left panels of fig. [5J There are 
no significant changes on the other cosmological parameters. 

This situation is also reproduced when the fastroll stage is included, that is when the fast roll transfer function 
D-jz(k) eq. (|3.35[) and fig. [3] is used, treating the scale k\ in eq. (|3.36|) as a MCMC parameter. 

That is, in the MCMC analysis with fast-roll included, we use the initial power spectrum eq. (|3.31|) which is modified 
by the fast-roll transfer function D-jiik). We computed once and forever D-jz(k) from eq. (|3.35[) [see fig. [3J. D-n{k) is 
a function of k and k\ with the scaling form eq. (|3.36p . ^>{x) being an universal function. 
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FIG. 7: Marginalized distributions (solid blue lines) and mean likelihoods (red dotted lines) for the parameters of the 
ACDM+sharpcut model. 



ACDMz | Cbni +sharpcut 


fa (best fit) 


r (best fit) 


r (95% CL) 


WMAP-3 only 


0.275 Gpc" 1 


0.150 


> 0.023 


CMB+LSS 


0.268 Gpc" 1 


0.051 


> 0.018 



TABLE II: The most likely values of fa and r and the lower bound on r in the ACDMz|Cbni + sharpcut) model. 



We then find 



minx|(ACDMz|C B Ni + fastroll) - minx|(ACDMz|C B Ni) 




(WMAP-3 only) 
(CMB+LSS) 



Correspondingly, in Table III we report the most likely values (ML) of k\ and r (we report also the best fit for the 
quartic coupling y for future use) where we used the marginalized probability in the (r, fa) plane as shown in the 
two right panels of fig. [5] Here ACDMz|Cbni + fastroll denotes the ACDMz|Cbni model with the fast roll Dn(k) 
included. 

We see a clear peak in y when fast-roll or a sharp cut are introduced in the CMB+SDSS fits. 

Wc sec that the gains in likelihood are more significant in the fast-roll case than in the sharpcut case. Clearly, this 
fit improvement through power modification by fastroll over power reduction by sharpcut is too small to constitue 
a real experimental evidence. But still, it is very interesting that the theoretically well founded approach based on 
fastroll works better than the purely phenomenological cutoff. This may be appreciated also from fig. [9] where the 
best fit for the Cj T multipoles are compared to the experimental data at low £. We see that the oscillatory form of 
the fastroll transfer function D-n(k), by depressing as well as enhancing the primordial power spectrum at long 
wavelengths, leads also to new superimposed oscillatory corrections on the multipoles. As far as fitting to current 
data is concerned, such corrections are more effective than the pure reduction caused by a sharp cutoff. 



ACDMz | Cbni +fastroll 


ki (best fit) 


r (best fit) 


r (95% CL) 


y (best fit) 


WMAP-3 only 


0.249 Gpc" 1 


0.146 


> 0.018 


0.031 


CMB+LSS 


0.266 Gpc" 1 


0.058 


> 0.018 


1.322 



TABLE III: The most likely values of fa, r and the quartic coupling y and the lower bound on r in the ACDMz | Cbni + fastroll) 
model. 
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We did not display in figs. [9] and [10] the ACDM+sharpcut results since they are indistinguishable from the 
BNI+sharpcut values. 

0.5 1 I 1 



0.4 




FIG. 8: Marginalized pair distributions in the (r, ki) plane, at 20%, 41%, 68% and 95% CL, with the WMAP-3 data alone 
in the top panels and with the CMB+SDSS data in the bottom panels. The left panels refer to the ACDMz|CBNi+sharpcut 
model, while the right ones are for the ACDMz|C*BNi+fastroll model, ki is in Mpc -1 . 

We plot in fig. [10] the best fit for the Cj E multipoles compared to the experimental data at low I. We see that 
for I — 2 and £ — 3 fast-roll and sharpcut models yield rather similar results (and better than the ACDM+r model) 
while for £ — 4 fast-roll produces a value closer to WMAP-3 data than sharpcut. 

We plot in fig. [TT] the Cf E multipoles computed in the best fit point to the WMAP-5 data compared to the 
experimental WMAP-5 data at low £. We see that both fast-roll and sharpcut models produce a reduction of the low 
EE multipoles including the EE quadrupole. 

Our fits imposing Cbni predict non-zero lower bounds on r: at 95% CL, we find that r > 0.023 when no cutoff is 
introduced, while r > 0.018 when either the sharp cutoff or the fastroll D(k) are introduced. The best fit values of 
the other cosmological parameters remain practically unchanged as compared to ACDM. Similarly their marginalized 
probability distributions are almost unchanged, with the natural exception of n s , which in BNI has a theoretical upper 
limit [see eq. (|5 ■ - 

C. Real Space Two Point TT-Correlator 

We display in fig. [12] the real space two point TT-correlation function C TT (9) for ACDM, sharpcut and fast-roll 
models, 

CTT W-^E(^ + 1 )^ TT fl(cos0). 

1=2 

We see that the ACDM correlator becomes really different from the two others only for large angles 8 > 1. Since 
all /-modes besides the lowest ones are practically identical in the three cases, this shows how dominant are the low 
multipoles in the large angle correlations. We also show the WMAP data, the width of the data is mostly due to the 
cosmic variance. 

As is clear from fig. rT2] both fast-roll and sharpcut models reproduce the two point correlator C TT (9) better than 
the pure slow-roll ACDM+r model. 
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FIG. 9: Comparison, with the experimental WMAP-3 data, of the theoretical Cj T multipoles computed in the best fit point 
of the various models of the main text. The error bars in the plotted range of £ are mostly due to cosmic variance. The insert 
contains an enlargement in linear scale of the first seven multipoles. BNI stands for binomial new inflation. The Cj T units 
are \p, K 2 ] and they are plotted as functions of the natural logarithm of £. Error bars of the WMAP-3 data are one-sigma 
(68%c.L). 



V. THE TOTAL NUMBER OF E-FOLDS OF INFLATION N to t = N + 6 ~ 66 

Another interesting observation is possible concerning the number N of efolds since horizon exit till the end of 
inflation. First of all let us clarify why in all our MCMC runs we keep N fixed. The reason is that the main physics 
that determines the value of N is not contained in the available data but involves the reheating era. Therefore, 
although technically possible, it is not reliable to fit N solely with the CMB and LSS data within a pure, near scale- 
invariant slow roll scenario. On the other hand, the quadruple depression allows to set an absolute wavelength scale 
for the primordial power, so we can check the consistency of our assumptions about N which fixes the total number 
of efolds of inflation. 

In the case of the ACDMz|CBNi+fastroll model with the CMB+LSS datasets, the most likely value of the quartic 
coupling y is slightly larger than unity. Then from fig. [3] we read a value ~ 14 for the ratio k\jm at horizon exit. 
It is important to compare the quadrupole mode scale fcg with the scale k\ that characterizes the fast-roll stage. 
The physical quadrupole (I — 2) wavemodes today &q are related to the particle horizon today 770 by 

k Q 770 = 3.342 

where the spherical Bessel function j 2 (k rjo) takes its maximun value, and 770 is given by Q 

3.29 



when one takes into account the acceleration of the universe expansion for < z < 2. Therefore, using the present 
value H [l8[ we obtain, 

k Q = 0.242 (Ope)" 1 . 
Notice that the value of kq is smaller than the characteristic scale k\ . 
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FIG. 10: Comparison, with the experimental WMAP-3 data, ol the theoretical C\ multipoles computed in the best fit point 
of the ACDM+r model, fast-roll and sharpcut models. Notice that for C™ and Cj E fast-roll and sharpcut models yield 
rather similar results (and better than the ACDM+r model), while for £ = 4 fast-roll produces a value closer to WMAP-3 than 
sharpcut. The Cj units are [fx K 2 ] and they are plotted as functions of the natural logarithm of i. Error bars of the WMAP-3 
data are one-sigma (68%c.L). 



k 


In a at horizon exit 


t v at horizon exit 


k Q = 0.242 Gpc" 1 


1.01 


0.0276 > 1/N 


fci = 0.266 Gpc" 1 


1.107 


0.0188 ~ 1/N 


k = 2 Gpc- 1 (WMAP) 


3.135 


< 1/N 


feo = 50 Gpc" 1 (CosmoMC) 


6.363 


< 1/N 



TABLE IV: The number of efolds since the beginning of inflation when the wavenumbers kQ, ki, ko exit the horizon. The 
quadrupole modes exit the horizon during the fast roll stage, about 1/10 of an efold before fast roll ends. k\ precisely exits 
the horizon at the transition from the fast roll to the slow roll stage. 



We display in Table IV the relevant wavenumbers: kQ, ki, kg and the number of efolds since the beginning of 
inflation when they exit the horizon. We see that the quadrupole modes exit the horizon during the fast-roll stage, 
approximately 1/10 of an efold before the end of fast roll. The mode k\ exit the horizon by lna = 1.107, very close to 
the point lna = 1.091 where e v — 1/N. That is, k\ precisely exits the horizon when fast roll ends and becomes 
slow roll. 

We denote by ko in Table IV the pivot wavenumbers in the WMAP H and CosmoMC codes [l3| , where the indices 
n s , r and the running of n s are computed. Both fco's exit the horizon well inside the slow roll regime. 

We read from Table IV that the total number of efolds of inflation is given by 

Ntot ~ N + 6 

since we have six efolds before the pivot wavenumber in CosmoMC exit the horizon followed by N efolds of inflation. 

We can compute the redshift 1 + zt since the begining of inflation till today comparing kq = 0.242 (Gpc) -1 (today) 
with k^ ltial = 0.910 k\ mtial = 12.7 m (at the begining of inflation). [Recall that 1 GeV = 1.564 x 10 41 (Gpc)" 1 ]. 
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FIG. 11: Comparison, with the experimental WMAP-5 data, of the theoretical Cf E multipoles computed in the best fit point 
to the WMAP-5 data of the ACDM+r model, fast-roll and sharpcut models as functions of £. Error bars of the WMAP-5 
data are one-sigma (68%c.L). Notice that both fast-roll and sharpcut models produce a depression of the low EE multipoles 
including the EE quadrupole. The Cf E units are \p K 2 ]. 



We use for m the value obtained from the scale of inflation [if 

M 2 



m 



M P 



where M is fixed by the amplitude of the scalar adiabatic fluctuations [1, [l5| . We obtain 

M = 0.57 x 10 16 GcV and m = 1.34 x 10 13 GeV for y = 1.322 . 
[Notice that these results are in agreement with the generic estimates eq. (|2.4j) .] Therefore, 

1 + 2% = 1.10 X 10 56 ~e 129 . 

Assuming a sharp transition from inflation to radiation dominated expansion, the redshift 1 + z\, can be written as 

1 



10 



-56 



l + Z b 



-N tot 



(5.1) 



where a r is the scale factor at the begining of the radiation dominated era and N to t is the total number of efolds 
during inflation (during fast-roll plus during slow-roll). 

The scale factor at the beginning (a r ) and the end (a eq , equilibration) of the radiation dominated era are related 

by 




where H and H eq stand for the Hubble parameter at the beginning and at the end of the radiation dominated era, 
respectively. For simplicity we assume instantaneous reheating in these formulas. 
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FIG. 12: The real space two point TT correlation function C TT (6) for ACDM, sharpcut and fast-roll models vs. the angle 9. 
The ACDM correlator differs from the two others only for large angles 6 > 1. Since all Z-modes besides the lowest ones are 
practically identical in the three cases, this shows how important are the low multipoles in the large angle correlations. Also 
shown are the WMAP data. The truly observed correlator runs approximately in the middle of the red band. The width of the 
data band is mostly due to the cosmic variance. The WMAP C TT (9) plotted here may not coincide, especially for the largest 
values 6 ~ n, with the correlator directly measured from sky maps due to the pixel weighting in the WMAP data analysis. 
The C TT {8) units are [/i K 2 ]. 



Furthermore, H eq and Hq are related by [if 



H e q — y2 Q m a e q 2 Hq . 

where Hq stands for the Hubble parameter today and tt m for the matter fraction of the energy density of the universe 
today. 

Using the current values of the cosmological parameters [IH we find 



1n -29 10- 4 MPI -67 10- 4 MPI (K0 , 

a r ~10 y d . (5.2) 



Inserting eq. (|5.2[) into ea. (|5.ip yields 

10" GeV " 6 ■ [b - 6) 

Since H must be below its value at the beginning of inflation ~ 10 14 GeV [see ea. (|2.5p ]. we conclude that 

N > 56 . (5.4) 

On the other hand, we know from BBN (Big Bang Nucleosynthesis) that H is at least larger than 1 MeV. This 
together with eq. (|5.3[) yields the upper bound 

N < 76 . (5.5) 

Furthermore, our MCMC simulations give good fits for N ~ 50 — 60. The bound ea. (|5.4p therefore favours N ~ 60 
which implies N to t ~ 66 and H ~ 3 x 10 10 GeV. In addition, from eqs. (|5.4p and (|5.5p we obtain the bounds 62 < 
Ntot < 82. 
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In summary, the fast-roll stage explains the quadrupole suppression and fixes the total number of efolds of 
inflation [I, [J. 

Our present MCMC analysis yields N to t ~ 66. More generally, the upper bound eq. (|5.5|1 implies N to t < 82. 

Changing N from 50 to 60 does not affect significatively the MCMC fits we present in this paper. This is partially 
due to the fact that a change on y can partially compensate a change on TV. Another hint to increase TV above 50 
comes from WMAP-5 that gives a larger n s and using the theoretical upper limit for n s [3 EH : 

1.9236... 

' N ' ^ ' 

which gives n s < 0.9679 . . . for N = 60. This value is compatible with the n s value from WMAP5+BAO+SN and no 
running 
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VI. APPENDIX 

As already shown in sec. IIV| the four types of correlation multipoles, TT, TE, EE and BB, among the CMB 
anisotropics can be numerically computed with good accuracy within the CosmoMC program starting from a given 
cosmological model with a fixed value for all its parameters. To be precise, this computation is performed by the CAMB 
subprogram, which is an evolution of CMB FAST [171 ]. CAMB can also compute, with several levels of approximations, 
the matter power spectrum observable today given the primordial power spectrum of density perturbations. 

On modern workstation CPU's, the calculation of two thousand scalar multipoles (related to the primordial curva- 
ture fluctuations) and one thousand tensor multipoles (related to primordial gravitational waves) takes less than one 
second. Thus, within a given type of cosmological model, the theoretical predictions for different choices of parameters 
can be produced at a very high rate. 

To test these predictions against the experimental data CosmoMC makes use of likelihood functions to assign 
different weights to different sets of correlation multipoles and matter power data. The experimental data and 
associated numerical code to evaluate such likelihoods is part of CosmoMC in the case of small-scale CMB experiments 
and LSS surveys. The data and likelihood routines for WMAP are not part of CosmoMC, but they can be downloaded 
from http: / /lambd aTgsfc.nasa.gov/| and integrated quite easily into CosmoMC, since its interface to CMB likelihoods 
has been appositely designed for WMAP. 

The WMAP likelihood code is particularly complex with respect to the other experiments, due to the wealth 
of experimental data, the variety of source of systematical errors and the importance of cosmic variance on lower 
multipoles. Indeed, it has significally changed and improved along the three WMAP releases. At any rate, it is used 
in CosmoMC exactly as released by the WMAP team. 

Finally, CosmoMC provides the MCMC engine, that is, the routines to perform suitable random walks (the chains) 
in the parameter space of a give cosmological model in such a way to reconstructs the experimental probability for 
the parameters from the distribution of values produced along the chains. In our simulations we always employed the 
default Metropolis rule, where the one-step transition probability from one set A of parameter values to the next is 
given by 

where g{\' , A) is a Gaussian proposal, or jump, probability and £(A) is the complete posterior likelihood, that is the 
product of the prior probability of choosing the starting point of the chain, times the likelihood obtained by comparing 
the theoretical prediction on multipoles (and matter power if LSS constraints are required) with the experimental 
data. As is well known from the theory of stationary Markov chains, in the limit of infinitely long chains no dependence 
is retained on <?(A', A) and the reconstructed profile is that of L(X) only. Of course, since actual chains have a finite 
length, suitable convergence tests are needed to verify that such a reconstruction is accurate enough. On pc-cluster 
running several parallel chains at once, CosmoMC offers very effective tests of this kind. 
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In order to test our own models based on Binomial New Inflation, with or without sharp-cut or fast-roll, we needed 
to modify some routines in CosmoMC. Since our changes with respect to the ACDM model are restricted to the 
primordial power spectrum, only routines relative to the so-called fast variables needed suitable modifications. Fast 
variables in CosmoMC are cosmological parameters that affect only the primordial spectrum so that, in a Monte Carlo 
step that proposes changes restricted to them, no need arises to perform the time-consuming recomputation of the 
transfer functions from a given primordial spectrum to observable multipolcs. Slow variables such as the baryonic and 
dark matter fractions, w^, w c , the optical depth r, and the present Hubble parameter Hq have the opposite definition. 
No change was done in the by far major portion of the CosmoMC program that deals with slow variables. 

We recall also that we let only the four slow variables mentioned above vary in our MCMC run (to be precise, we 
used the default choice of CosmoMC which replaces H with 9, the ratio of the approximate sound horizon to the 
angular diameter distance), keeping all other slow cosmological parameters, such as the parameter of the dark energy 
equation of state or the neutrino density fraction, to the values of the standard ACDM model. 
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